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I. INTRODUCTION 



CP violating asymmetry is one of the most important areas in the decays of bot- 
tom hadrons. In the standard model(SM), a non-zero complex phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix is responsible for CP violating phenomena. In recent 
years CP violation in several B decays such as B° — > J/ipK^ and B° — > K + ir~ has indeed 
been found in experiments [j], 0]. Due to its much higher statistics, the Large Hadron 
Collider (LHC) will provide a new opportunity to search for more CP violation signals. 

Direct CP violating asymmetries in 6-hadron decays occur through the interference of 
at least two amplitudes with the weak phase difference and the strong phase difference 
5. The weak phase difference is determined by the CKM matrix while the strong phase is 
usually difficult to control. In order to have a large CP violating asymmetries signal, we 
have to apply some phenomenological mechanism to obtain a large 5. It has been shown 
that the charge symmetry violating mixing between p° and u can be used to obtain 
a large strong phase difference which is required for large CP violating asymmetries. 
Furthermore, it has been shown that the measurement of the CP violating asymmetries 
can be used to remove the mod(7r) ambiguity in the determination of the CP violating 
phase angle aj3|-0|- 

Naive factorization approximation has been shown to be the leading order result in the 
framework of QCD factorization when the radiative QCD corrections of order 0(a s (rrib)) 
(rrif, is the 6-quark mass) and the 0(l/mb) corrections in the heavy quark effective theory 
are neglected jsj. In naive factorization scheme, the hadronic matrix elements of four- 
quark operators are assumed to be saturated by vacuum intermediate states. Since the 
bottom hadrons are very heavy, their hadronic decays are energetic. Hence the quark pair 
generated by one current in the weak Hamiltonian moves very fast away from the weak 
interaction point. Therefore, by the time this quark pair hadronizes into a meson, it is 
already far away from other quarks and is unlikely to interact with the remaining quarks. 
This quark pair is factorized out and generates a meson [9], [Ioj|. This approximation 
can only estimate the CP violation order neglecting QCD correction. Furthermore, as 
pointed out in previous studies 0-0|, in order to taken into account the nonfactorizable 
contributions, an effective parameter, N c , is introduced. The deviation of the value of N c 
from the color number, 3, measures the nonfactorizable effects in the naive factorization 
scheme. Obviously, N c should depend on the hadronization dynamics of different decay 
channels. In this scheme, CP violation depends strongly on N c values, which makes the 
results uncertainties. 

In the heavy quark limit, QCD factorization jsj includs nonfactorization strong interac- 



tion correction, and the decay amplitudes can be calculated at leading power in and 
at next-to-leading order in a s , which can be expressed in terms of form factors and meson 
light-cone distribution amplitudes. One can take into account the nonfactorizable and 
chirally enhanced hard-scattering spectator and annihilation contributions which appear 
at order 0{a s {mi ) )) and 0(l/mb), respectively. In this work we adopt the QCD factor- 
ization scheme including order-a s correction to compute CP violating asymmetry of the 
decay B° s — > K°7t + 7t~ via the p — u mixing mechanism. As will be shown later, the CP 
violating asymmetries in this decay channel could be large and may be observed in the 
LHC experiments. 

The remainder of this paper is organized as follows. In Sec. II, we present the form 
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of the effective Hamiltonian and the general form of QCD factorization. In Sec. Ill, we 
give the formalism for CP violating asymmetries in B° s — > K°7t + tt~ decay. In Sec. IV, we 
calculate the branching ratio for decay process of B° s — > K°p°(uj) via p — u mixing. We 
briefly discuss the input parameters in Sec. V. The numerical results are given in Sec. VI. 
In Sec. VII we discuss the possibility to observe the predicted CP violating asymmetries 
at the LHC Summary and conclusions are included in Sec. VIII. 



II. THE EFFECTIVE HAMILTONIAN 



With the operator product expansion the effective Hamiltonian in bottom hadron 
decays is 



v p=u,cq=d,s 
10 

+ ^ C i°i + C 77°7 7 + C 8g 8g ] + H.C, (1) 
i=3 

where Cj {i = 1, 10, 7j, Sg) are the Wilson coefficients, V p b, V pq are the CKM matrix 
elements. The operators 0« have the following form: 

0{ = F7^(l - 75)&n M (l - 7s)p, O p 2 = ^7^(1 - 7 5 )6/3g /3 7 / "(l - 7s)p«, 

3 = q lfl (l - 7s)& £ </7 M (l - 75)g', 4 = q a ^(l - 7s)bp £ ? p -f(l - 7s)^, 

q' q' 

5 = g T „(l - l5 )b £ q>r{l + J5)q', 6 = q al ,{l - 75)6,3 £ q'^l + l5 )q' a , 

q' q' 

7 = - 75)& £ e q/ q'r(l + JsW, 8 = §g a7/i (l - 75)6/3 £ e^T^l + 75)^, 

9' g' 

9 = |?7m( 1 -75)^£e (? 'g , 7 At (l -7s)?', O w = ffeT^ 1 ~ 7s)&/3 £ V^aT^l - 75)^, 

<?' 9' 
77 = + 75)^6, 8g = =^m b so^{\ + 75)^6, 

(2) 

where a and /3 are color indices, 0\ and Of are the tree operators, O3 — 0@ are QCD 
penguin operators which are isosinglets, O7 — 0\q arise from electroweak penguin operators 
which have both isospin and 1 components. and Og g are the electromagnetic 

and chromomagnetic dipole operators. e q i are the electric charges of the quarks and 
q' = u, d, s, c, b is implied. 

The Wilson coefficients can be calculated at a high scale My/ and then evolved to scale 
rrib using renormalization group equation. In QCD factorization, We consider weak decay 
B s — > M\M.2 (Mi, M2 refer to K° and p° mesons, respectively) in the heavy-quark limit. 
Up to power corrections of order AQco/^b, the transition matrix element of an operator 



3 



Oi in the weak effective Hamiltonian is given by[8j] 

(M 1 M 2 \O i \B) = Y^F^ Ah {m 2 2 ) f duTt(u)$ Ma (u) 

+ f dZdudvT l II (^U,v)$ B (0®M 1 (v)®M 2 (u) 

Jo 

if Mi and M 2 are both light, (3) 

Here 1,2 (m^) denotes aB-> M 12 form factor, and is the light-cone distri- 

bution amplitude for the quark-antiquark Fock state of meson X. T^{u) and T/ 7 (£, u, u) 
are hard-scattering functions, which are perturbatively calculable. The hard-scattering 
kernels and light-cone distribution amplitudes (LCDA) depend on a factorization scale 
and scheme, which is suppressed in the notation of ([3]). Finally, mi |2 denote the light 
meson masses. 

We match the effective weak Hamiltonian onto a transition operator, the matrix ele- 
ment is given by (X P D ^ = V p bV* D with D = d or s) 

(M[M' 2 \U c z\B) = ^ (KK\U +1I\B). (4) 

p=u,c 

Using the unitarity relation 

XL D) + >S» + Ai D) = (5) 

we can get 



E T'tx = E 



p=u,c 



p=u,c 



S pu a 1 (M 1 M 2 ) A{[q s u][uD\) + 5 pu a 2 (M 1 M 2 ) A([q s D][uu}) 



+ W 



K(M!M 2 ) - aKM.M,)) A([q s q][qD]) + (a« EW (MiM 



V - «4,EW 



(MiM 2 )) ^£-e q A([q s q)[qL 



- X 



+ a 



a 



;(MiM 2 ) Y,A([q 8 D][qq])+al(M 1 M 2 ) ^ A([g.g][gZ>]) + ^(M^) ^ - e 9 



4,EW 



(MiM 2 )X;^e ? A([^][gL>]) 



where the sums extend over q = u,d, s, and g s denotes the spectator antiquark. The 
operators A([q Ml q Ml ][q M2 q M2 \) also contain an implicit sum over q s = u,d,s to cover all 
possible 5-meson initial states. 

Next we need change the annihilation part 
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p=u,c 



E = E ^ 

p=u,c 

5 pu b 1 (M 1 M 2 ) Y^B([uq'][q'u][Db]) 

q' 

+ 8 pu b 2 (M l M 2 ) ^B([uj][ifD][ub])) 



A 



(£>) 



b 3 (M 1 M 2 ) Y,B{[qq'][q'D][qb}) 



q,g' 



+ 6 4 (M 1 M 2 ) 



+ 6 4| ew(M 1 M 2 ) ^-e, J B([gg'][g'g][ J D6]) 



(7) 



where fr^EW and B are given by following. The coefficients of the flavor operators 
can be expressed in terms of the coefficients af defined in [s| as follows: 

a x {M x M 2 ) = a 1 {M 1 M 2 ), 
a 2 (M 1 M 2 ) = a 2 (M 1 M 2 ) , 



al{M x M 2 ) 
a p iEW {M 1 M 2 ) 



^(MiM^+ogCMiMa); 
if MxM 2 = PV, 

if M!M 2 = PV , 

ag(MiM 2 ) + a?(MiM 2 ) ; 
if MxM 2 = Py, 

a?o(M 1 Af 2 )+rJ f »ag(M 1 M 2 ); 
if MxM 2 = P^ , 



For pseudoscalar (P) meson M 1; the ratios r^ 1 are defined as 



2m^ 



m & (/i) (m q + m s )(n) 



(9) 



All quark masses are running masses defined in the MS scheme, and m g denotes the 
average of the up and down quark masses. For vector (V) meson M 2 we have 



(10) 
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where the scale-dependent transverse decay constant fy is defined as 

(V(p,e*)\qa ia ,q'\0) = f^et - Pv e*) . 



Ill 



Note that all the terms proportional to r^ 2 are formally suppressed by one power of 
Aqcd l m b in the heavy-quark limit. 

The general form of the coefficients af at next-to-leading order in a s is 



a\{M x M 2 ) 



+ N c 4vr 



4vr 2 

V i (M 2 ) + —H i (M 1 M 2 ) 



(12) 



where N c is the number of colors, the upper (lower) signs apply when % is odd (even). It 
is understood that the superscript 'p' is to be omitted for i = 1,2. The quantities Vi(M 2 ) 
account for one-loop vertex corrections, Hi{M\M 2 ) for hard spectator interactions, and 
P[(MiM 2 ) for penguin contractions. The Ni(M 2 ) and Cf are given by 



N t {M 2 




1 



C7 



% = 6, 8 and M 2 = V, 
all other cases. 

iV?-l 



The vertex corrections are given by 



Vi(M 2 



Jg 1 dx$ M2 (x) 
(i = 1-4,9, 



2iV; 



12 In ^ - 18 + ^(x) 



0), 



jl dx §m 2 {x) y - 12 In St + 6 - #(1 - s ) 
(i = 5,7), 



Jo 1 rfX 

(i = 6,8), 



6 + /i(a;) 



with 



,'1-20;, 
3 ( — In x — m ) + 



+ 



1 — x 
2 In a; 



1 — x 



(3 + 2ivr)lnx 



2Li 2 (a;) - ln 2 x 
(x ■<->■ 1 — x) 



(13) 
(14) 



(15) 



(16) 



h(x) = 2 hi 2 (x) — \n 2 x — (1 + 27ri) In x — (x •<-)■ 1 — x) . 



;i7) 



The constants —18, 6, —6 are scheme dependent and correspond to using the NDR scheme 
for 75. The light-cone distribution amplitude (LCDA) $m 2 is the leading-twist amplitude 
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of M 2 , whereas $ m2 is the twist-3 amplitude. LCDA for pseudoscalar and vector mesons 
of twist- 2 are 



$p(x,ii) = 6x(l — x) 
^J(x,/jl) = 6x(l — x) 
$I(x,/-0 = 6x(l-x) 



n=l 

oo 

n=l 

oo 

i + £^ y (^/ 2 (2*-i) 



n=l 



(18) 



and twist-3 ones 



%(x) = 1, $ a (x) = 6x(l-x), 



$ v (x,fx) = 3 



n=l 



(19) 



where C n {x) and P n (x) are the Gegenbauer and Legendre polynomials, respectively. a n (n) 
are Gegenbauer moments that depend on the scale \x. §^_(x,fi) and ^J(x,/i) are the 
transverse and longitudinal quark distributions of the polarized mesons. 

At order a s a correction from penguin contractions is present only for i = 4, 6. For 
i = 4 we obtain 



4nN, 

+ c 3 



4, m 6 2 

3ln- + --G M2 ( Sp ) 



|ln^ + |-G Wa (0)-G Ma (l) 



+ (C 4 + C 6 



^ln^-( n/ -2)G Ma (0) 



— Gm 2 (s c ) — Gm 2 (1) 



2r;:; /' ; -^-<i-u,(.r); 



o 1 - x 



(20) 



where n/ = 5 is the number of light quark flavors, and s u = 0, s c = (m c /m b ) 2 are mass 
ratios involved in the evaluation of the penguin diagrams. The function G M2 (s) is given 
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by 



G M2 {s) 



G(s, x) 



/ dx G(s — ie,l — x) <&m 2 (x) 
Jo 



—4 / du u(l — u) ln[s — u(l — u)x] 
Jo 

2(12s + 5x-3xliis) 



(21) 



9x 



Ay/As - x (2s + x) 



arctan 



x 



3x 3 / 2 \/ 4s — x 

For % — 6, if M 2 is a vector meson, the result for the penguin contribution is 

C F a s 



(22) 



AnN r 



C\ Gm 2 ( s p) + G 3 Gm 2 (0) + G ! m 2 (1 



+ (c 4 + c 6 



(n f -2)G M2 (0)+G M2 (sc 



(23) 



In analogy with ( I2T]) . the function G M2 (s) is defined as 



Gm 2 ( s ) 



dx G(s — ie, 1 — a;) $ m2 (x) 



(24) 



Electromagnetic corrections are present for i — 8, 10 and correspond to the penguin 
diagrams. For i = 10 we obtain 



P?o(M 2 



a 



9nN r 



(d + iV c C 2 ) 



4 , m 6 
- In — 

3 ji 



+ g - Gm 2 {s p ) 



For i = 8 



Pf(M 2 



3C5f Z 1 ^*^) 



(Ci + iV c C 2 )G M2 (s p ), 



(25) 
(26) 



9vriV c 

if M 2 is a vector meson. 

The correction from hard gluon exchange between M 2 and the spectator quark is given 

by 



H i (M 1 M 2 ) 



1 /•! 



dx / dy 



B MlM2 rriB 
Am 1 m 2 *b Jo -hi 

Aii 



(27) 
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for i = 1-4,9,10. 



H i (M 1 M 2 ) 



BmiM 2 m B 
Am 1 M 2 Jo 



dx I dy 



§m 2 {x)<&mM 
xy 



_l_ r 



$ M2 (x)$ mi (y) 
xy 



(2f 



(29) 



for i = 5, 7, and H i (M 1 M 2 ) = for z = 6, 8. 
where A# is defined by 

Jo 4 

with $b(£) is one of the two light-cone distribution amplitudes of the B meson. 

If Mi = P, M 2 = V, f refers to decay constant of relevant meson, A MlMi and B MlM2 
are given by 

A MlM2 =i^(-2)m Ml e* M -p B F*-+ Ml (0)f M2 , (30) 



v/2 



G F 

Bm!M 2 — — — ^/b s /mi/m 2 - 



(31) 



where and are the mass and polarization vector of the vector meson. F ~^ 1 is 
the form factor for B — > Mi transition. 

We recall that the term involving r^f 1 is suppressed by a factor of Aqcd/^6 in heavy- 
quark power counting. Since the twist-3 distribution amplitude $ mi (y) does not vanish 
at y — 1, the power-suppressed term is divergent. We extract this divergence by defining 



a parameter X H 1 through 

"dy 



Jo y Jo y 



$ mi (y) - $ mi (l) 
1 



(32) 



The remaining integral is finite (it vanishes for pseudoscalar mesons since $ p (y) = 1), 
but is an unknown parameter representing a soft-gluon interaction with the spec- 
tator quark. Since X^ 1 varies within a certain range (specified later) and Xjf ~ 
ln(mfe/AQCD)|8(, we treat the resulting variation of the coefficients as an uncertainty. 
We also assume that X H X is universal, i.e., that it does not depend on Mi and onjthe 
index i of Hi(MiM 2 ) . For the convolution integrals, one can find the results in Ref. 
For the annihilation contribution, one can get[sj]: 



C 3 Ai + C 5 (Ai + Ai)+N c C 6 A f 3 



(33) 
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°3,ew - m 



C 9 A\ + C 7 (Ai + A f 3 ) + N C C 8 A( 



(34) 



The weak annihilation kernels exhibit endpoint divergences, which we treat in the 
same manner as the power corrections to the hard spectator scattering. The divergent 
subtractions are interpreted as 



1 dy ^ Y^i 

o y 



Iny 



dy—^-~{X^f 

o y z 



(35) 



and similarly for M 2 with y — > x. The treatment of weak annihilation is mo del- dependent 
in the QCD factorization approach. We treat X A as an unknown complex number of order 
ln(mft/AQCD) and make the simplifying assumption that this number is independent of 
the identity of the meson Mi and the weak decay vertex. Here, 



+ rfrf (X 2 A -2X A ) 



3 X A -4 + 



7T 



(36) 



A\ * 



Qua, 



— 3r 



+ 4) +rf (XI 



Ai 2 



X A — 2X A 



2X A + — 



71' 



(37) 



y4g w — 67TQ; S 



3rf 2 (2X A -l)(2-X A ) 



- rf (2X1 " *a) 



(3? 



and A{ = A{ = 0. Here, Mi is fT° meson and M 2 is p° meson. 



III. CP VIOLATION IN B° s K°7r+vr DECAY 
A. Formalism 

In the vector meson dominance model [12j, the photon propagator is dressed by cou- 
pling to vector mesons. Based on the same mechanism, p — u mixing was proposed (l3| . 
The formalism for CP violation in the decay of a bottom hadron, B s , will be reviewed in 
the following. The amplitude for B s — > K°7r + 7T~ , A, can be written as 

A = (tt + tt-K°\H t \B s ) + (tt + tt-K°\H p \B s ), (39) 

where H T and H p are the Hamiltonians for the tree and penguin operators, respectively. 
We define the relative magnitude and phases between these two contributions as follows: 

A = (n + n-K \H T \B s )[l + re i5 e% (40) 
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where 5 and are strong and weak phase differences, respectively. The weak phase 
difference arises from the appropriate combination of the CKM matrix elements: <fi = 
aig[(VtbV t *) I \V u bV* s )}. The parameter r is the absolute value of the ratio of tree and 
penguin amplitudes, 



(ti + 7i-K \H p \B s ) 



(tt+tt-K°\H t \B s ) 
The amplitude for B s — > K°7t + tt~ is 

A = (tt + tt-K°\H t \B s ) + (tt + tt-K°\H p \B s ). 
Then, the CP violating asymmetry, a, can be written as 

— 2r sin 5 sin <b 



\A 


2 _ 


A\ 2 


\A 


2 +|A| 2 



1 + 2r cos 5 cos 6 + r 2 



(41) 



(42) 



(43) 



We can see explicitly from Eq. (40) that both weak and strong phase differences are 
needed to produce CP violation, p — u mixing has the dual advantages that the strong 
phase difference is large and well known 0, 0] . In this scenario one has 



(n + n-K°\H T \B s ) 



(44) 



(ir+TT- K°\H P \B S ) = J^fL^+p^ + ^+p:) 



(45) 



where tv(V — p or uj) is the tree amplitude and py is the penguin amplitude for producing 
a vector meson, V. ty(V = p or uj) is the tree annihilation amplitude and p v is the penguin 
annihilation amplitude. g p is the coupling for p° — > 7r + 7r~, U puJ is the effective p—u mixing 
amplitude, and sy is from the inverse propagator of the vector meson V, 

Sy = s - m v + imyVy, (46) 

with yfs being the invariant mass of the 7r + 7r~ pair. 

The direct uj — > 7r + 7r~ is effectively absorbed into H puj , leading to the explicit s de- 



pendence of H puJ [14J. Making the expansion U puj (s) = il pa; (rrt 2 J ) + (s — rn^)W 



the p — u mixing parameters were determined in the fit of Gardner and O'Connell |15| : 
Retl pul {ml) = -3500 ± 300 MeV 2 , Imfi pt ,(m 2 ) = -300 ± 300 MeV 2 , and n^(m 2 ) = 
0.03 ± 0.04. In practice, the effect of the derivative term is negligible. From Eqs. 
(40) (41) (44) (45) (46) one has 

re i6 e i4> = npo,(^+K) + ^(Pp+P°) 
ILpuitv + t*) + S u (t p + tf) ' 



11 



Defining 

P^±P^ = r > e i(6 v U) hi±% = ae is Pp±A = 0e iS, (48) 
tp + t* P t P + t a P P«+P£ 

where S a , 5p, and S q are strong phases, one finds the following expression from Eq. (47): 

re = r e "— — J — . (49) 

ae lSa , /3e i<5/3 , and re lS will be calculated in the QCD factorization approach later. With 
Eq. (49), we can obtain rsin5 and rcos<5. In order to get the CP violating asymmetry, 
a, in Eq. (43), sin0 and cos0 are needed. is determined by the CKM matrix elements. 
In the Wolfenstein parametrization [l6j], one has 

sin0 = , (50) 

y/\p(l ~ p) -7] 2 ] 2 + T] 2 

cos = . (51) 

y/\p(l -p)~ V 2 ] 2 + V 2 



B. CP violation via p — uj mixing 

In the following we will study the CP violating asymmetries in the following decay: 
B° s -»■ K°p°(u) -»■ K°ir + ir~. With the Eq. (4)(6)(7)(8), we can calculate the decay am- 
plitudes in QCD factorization scheme. The expressions for the — > K°p°(u) amplitudes 
are given by 

V2A§o^ K o p o = A K o p o(5 pu a>2 — ol\ + ~a% EW 

+ \*Iew-PI + \pI,ew)i (52) 



V2A B o^ KOuJ = A K o OJ (d pu a 2 + 2o| + a\ + \-a\ 



where 



1 1 



A^opo = (-2)z%m p o£* • PB F o B ^ K °(0)f pO , (54) 



A^o, = (-2)i^=m u el ■ p B F B ^ K °(0)U (55) 
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Here Fq denote B s — > K° meson form factor. m p o, are the mass of p° and ou mesons. 
£* , £* correspond to polarizing vectors. / refers to the decay constant. Then we can get 



V2A 



B s ^K°p° 



iKOpO 



Sp U a 2 ,K0p0 - a p 4 K0p0 - 7 X al KOp0 



3 1 



(56) 



V2A 



A K °u 



(57) 



where the form of the coefficients af at next-to-leading order in a s is given by Eq.(12), 
which Mi is K° meson and M 2 is p° meson, fa is the weak annihilation contribution in 
QCD factorization. 7 X is chirally-enhanced terms which we have denoted above. 
From Eq. (6) (7) (48), one can get 



ae 



t„ + K 
t P + t° 



Qi 

Q2 



(58) 



+ 7 X l a 8,i^ - a 8,K°u>) 



(59) 
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Q 2 — A K o p o< Sp U a 2t K°p ~ ( a 4,K°p° ~ a 4,K»p") 



+ a 7,K°p° ~~ a 7,K»p») + 



a 10,i^°p° ~ a 10,i^V 



+ 7 X l a 8,i^V ~ a 8,K°p°, 



(60) 



In a similar way, with the aid of the Fierz identities, we can evaluate the penguin 
operator contributions p p and p w . From Eq. (48) we have 



(3e 



(61) 



where 



Q3 = A K o p o 



c kP c 

a A,K»pO ~ lx a Q,K0p0 



3 1 



+ 7x °8,X'Vo) ~ 03 + ~/?3,£W 



(62) 



and 



where 



Q 4 



~ 7;( a 10,K°u + 7^ag iX o w ) + /5 3 - ~f3 3 ,EW 



r ' e i(S q +4>) 



r ' e iS q 



tp+n ~ q 2 ' 



Qa 



v tb v t * d 



v ub v: d 



v tb v t * d 



V ubV* d 



a/[p(1 - p) - T] 2 ] 2 + r] 2 
(1- \ 2 /2)(p 2 + r] 2 ) 



(63) 

(64) 
(65) 

(66) 
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It can be seen that r' and S q depend on both the Wilson coefficients and the CKM 
matrix elements, as shown in Eqs. (65). Substituting Eqs. (58) (61) (65) into Eq. (49), 
we can obtain r, sin 8, and cos 5. Then, in combination with Eqs. (50) and (51) the CP 
violating asymmetries can be obtained. 



IV. BRANCHING RATIO OF -> K°p°(u) 



The matrix element for B s — > P and B a — > V (where P and V denote pseudoscalar 
and vector mesons, respectively) can be decomposed as follows [17j : 



(P\J„\B S ) = [p Bs+PP - 



m 



B a 



m 



k 2 



-k ) F^k 2 ) 



+ 



m% s - mp 



k 2 



KF {k 2 



{V\J,\B S ) = e^e w p p Ba P v V(k 2 ) 

m Bs + m v 

e*(m Bs + m v )A 1 (k 2 ) 6 — 

m Bs + m v 

e* ■ k 

x (Pb s +Pv)»Mk 2 ) ~ -j^- 2m v ■ k^A z (k 2 ^ 

+ i e -g^2m v ■ k^A (k 2 ), (67) 

where J M is the weak current (J M = g7 M (l — 75)6 with q = u,d,s), 
PB s ( m Bs),Pp( m p),Pv( m v) are the momenta (masses) of B S ,P,V, respectively, k = 
Pb s ~Pp{Pv) f° r B s — > P(V) transition and e M is the polarization vector of V. Fi {i = 0, 1) 
and Ai [i = 0, 1, 2, 3) in Eq. (67) are the weak form factors which satisfy i*i(0) = -^0(0), 
A 3 (0) = 4,(0), and A 3 (k 2 ) = [(m B + m v )/2m v }A 1 (k 2 ) - [(m B - m v )/2m v ]A 2 (k 2 ). 

With the factorizable decay amplitudes in Eq. (56) (57) we can calculate the decay 
rate for B s to a pseudoscalar meson (P) and a vector meson (V) transition by using the 
following expression [la ]: 

T(B S ^PV) = ^-r\A(B s ^PV)/(e- PBs )\ 2 , (68) 

07T illy 

where 



Pc 



\l\ m B s - ( m P + m v) 2 }[m 2 Bs - (m P - m v ) 2 } 



2m 



B s 



is the cm. momentum of the product particle and A(B S — > PV) is the decay amplitude. 
In the QCD factorization approach. Here V^ ,p are the CKM factors, 

Vl=\V ub V: q l for * = 1,2, (69) 
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and 



V u y = \V tb V t * q \, for i = 3, ....,10. (70) 

In our case we take into account the p — u mixing contribution when we calculate the 
branching ratios since we are working to the first order of isospin violation, we can 
explicitly express the branching ratio for the process B s — > K°p°(uj) as the following: 



BR{B S K°p°{u)) 

■|M(a 1) a 2 )-7Xo(a 3) ..,fl 10 )] 



p 

+ [V^Al(a u a 2 ) - V^A^{a 3 , a 10 )] 

x 7 " P l • r I 2 , (71) 

[s - mi) + imj: u 



where Y Ba is the total decay width of B s . 



V. INPUT PARAMETERS 



In the numerical calculations, we have several parameters, i.e. N c and the CKM matrix 
elements in the Wolfenstein parametrization. For the CKM matrix elements, which should 
be determined from experiments, we use the results of Ref. 



A 



B _ n i 09+O.022 

0.2253 ± 0.0007, 



f) = 0.341 ±0.013, 
A = 0.808^. 



(72) 



In QCD factorization scheme, since power corrections have been considered, N c is only 
color parameter, hence we use N c = 3. In naive factorization iV c includes the nonfatoriz- 
able effects which are model and process dependent and cannot be theoretically evaluated 
accurately and can be determined by experiment. 

The running quark masses is taken by the scale /1 in B s decay. One has 



m b (m b ) = 4.2GW, m c (m b ) = 0.91GeV, 
m u (m b ) = m d (m b ) = 0, m s (2AGeV) = 0.095GW. 

The values of the scale dependent quantities fyiu) and 0^2 (aO are given for p 
The value of Gegenbauer moments are taken from [19] . 



(73) 
lGeV. 



a 



a 



a 1 
1 



, A" 



fp 
f U 



0. 

0. 

0, 

at" = 0, 
= 0.06 ±0.03, 
216 ± 3MeV, 
187 ± 5MeV, 



a 



a 



a 



a. 



±U1 



0.15 ±0.07 
0.15 ±0.07 
= 0.14 ±0.06 
= 0.14 ±0.06 
a% = 0.25 ±0.15 
f^(p) = 165±9MeV, 
ti(p) = 151±9MeV, 



(74) 
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For B s meson, we use the value [2]: 

r = lA7ps, m Bs = 5.366GW (75) 



The Wilson coefficients q can be found in [8j . As discussed in detail in |8j , there are large 
theoretical uncertainties related to the modeling of power corrections corresponding to 
weak annihilation effects and the chirally-enhanced power corrections to hard spectator 
scattering. So we parameterize these effects in terms of the divergent integrals Xh (hard 
spectator scattering) and Xa (weak annihilation). We also model these quantities by 
using the parameterization [8|] 

X A = (l + ^ e ^) m !j£. A ft = .5GeV, (76) 

Aft 

and similarly for Xh- Here <£a is an arbitrary strong- interaction phase, which may be 
caused by soft rescattering. The fitted qa and if a are taken from [20|. For B s — > PV 
decay, p p A v ~ 0.87, $ P A V ~ —30°. For the estimate of theoretical uncertainties, we shall 
assign an error of ±0.1 to pa and ±20° to 6& 20l . 



The form factors associated with the weak transitions depend on the inner structure 
of the hadrons and are hence model dependent. Here we will consider the form factors 
obtained in several phenomenological models. For B s decay form factors, we will use the 
results (form factors are referred to the ones at q 2 = 0): 

1). Model 1 [8] 



pBs^K = o 32 ± Q.05, 



2). Model 2 (in the pQCD approach) [2l| 



^B s —>-K _ n 9 4 +0.05+0.00 
— u - z ^-0.04-0.01' 



3). Model 3 (form factors obtained by QCD sum rules) |22j| 



iBs^-K _ n qn+0.04 
— u - ou -0.03' 



4). Model 4 (light-cone sum rule calculation based on heavy quark effective theory) j23j 



F B.-nc = o.296 ± 0.018, 



5). Model 5__(A light cone quark model in conjunction with soft collinear effective 
theory) [24| 



F^ K = 0.290, 
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FIG. 1: Plot of a as a function of y/s corresponding to central parameter values of CKM matrix 
elements for B° s -> K°p°(ui) -»■ K°ir + ir~. 

6). Model 6 (lattice QCD calculation) |25| 



f b.-+k = g 23 ± 0.05 ± 0.04. 

In above Models, the k 2 dependence of the form factors has the following form under 
the nearest pole dominance assumption: 

h(k 2 ) = (77) 

where h could be F , and m h is the pole mass. 

It is noted that since the value of k 2 (which is actually the square of the mass of the 
factorized light meson) is much smaller than the square of the pole mass which is of order 
mf, only the values of the form factors at k 2 = are most relevant and hence how the 
form factors depend on k 2 has little effects (less than 2%). From the above values we see 
that the form factor B s — > K at q 2 = ranges from 0.23 to 0.31. 



VI. NUMERICAL RESULTS AND DISCUSSIONS 
A. CP violation via p — uj mixing in B® — > K°ir + ir~ 

In the numerical calculations, we find the CP violating asymmetry, a, is large when 
the invariant mass of tt + tt~ is in the vicinity of the u resonance within QCD factorization 
scheme. 

In the respective error ranges, when yfs = 0.782 GeV, we get maximum CP violating 
asymmetry 

a = (45.9tlf 7 ±^) x 10- 2 (78) 
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s" 2 (GeV) 

FIG. 2: Plot of sin 5 as a function of y/s corresponding to central parameter values of CKM 
matrix elements with p — uj mixing for B® — > K°p°(oj) — > K ir + ir~. 




073 0.74 0.75 0.76 0.77 0.78 0.79 80 0.81 0.82 



s (GeV) 

FIG. 3: Plot of r as a function of y/s corresponding to central parameter values of CKM matrix 
elements with p — u mixing for B® — >• K°p°(uj) — > K°tt + tt~ . 

In QCD factorization, the theoretical errors are large which follows to the uncertainties 
of results. Generally, power corrections beyond the heavy quark limit give the major the- 
oretical uncertainties. This implies the necessity of introducing l/m^ power corrections. 
Unfortunately, there are many possible l/m b power suppressed effects and they are gen- 
erally nonperturbative in nature and hence not calculable by the perturbative method. 
There are more uncertainties in this scheme. The first error refers to the variation of the 
CKM parameters. The second error comes from form factors and decay constants. The 
third error corresponds to the Gegenbauer moments. The last error is the wave function 
of the B s meson characterized by the parameter A#, the power corrections due to weak 
annihilation and hard spectator interactions described by the parameters pa,h, 4>a,h, re- 
spectively Using the central values of above parameters, we first calculate the numerical 
results of CP violation and branching ratio, and then add errors according to standard 
deviation. In Fig.l, We give the central value of CP violating asymmetry as a function of 
y/s. From the figure one can see the CP asymmetry parameter is dependent on y/s and 
changes rapidly due to p — u> mixing when the invariant mass of h + ti~ is in the vicinity of 
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the u resonance. The CP violating asymmetry vary from around —37% to around 45%. 

From Eq. (43), one can see that the CP violating asymmetry parameter depends on 
both sin S and r. The plots of sin S and r as a function of \fs are shown in Fig. 2 and Fig. 
3. It can be seen that when p — u mixing is taken into account sin 6 and r change sharply 
when the invariant mass of 7r + 7r~ is around 0.782 GeV. From the Fig. 2, one can find 
p — u mixing make the sin S value oscillate from —0.56 to 0.44 which can not reach the 
value —1. This result is not in agreement with conclusion from naive factorization which 
can measure the CP violating parameter to remove the mod(7r) phase uncertainty in the 
determination of the CKM angle a arising from the conventional determination through 
sin 2a Q. 

We have shown that p — u mixing does enhance the direct CP violating asymmetries 
and provide a mechanism for large CP violation in QCD factorization scheme. On the 
other hand, it is important to see whether it is possible to observe these large CP violating 
asymmetries in experiments. This depends on the branching ratio for B® — > K°p°(u). 
We will study this problem in the next section. 



B. Branching ratios via p — u mixing in B® — > K°p°(u) 

Including p — u mixing, we calculate the values of branching ratios for B° s — > K°p°(u). 
Base on the reasonable parameters range, we obtain the branching ratio of — > K°p°(u) 
is (9. 8^2-0 7) x 10~ 7 which is consistent with the result [lo[. In other words, although we 
calculate the branching ratio due to p — u mixing in QCD factorization scheme, we find 
the contribution of p — u mixing for branching ratio is small and can be neglected, p — u 
mixing mechanism presents new phase differences and produce extremely small effect for 
branching ratio of B° s — > K°p°(u). 

VII. DISCUSSIONS ON POSSIBILITY TO OBSERVE CP VIOLATING ASYM- 
METRIES AT THE LHC 

The LHC is a proton-proton collider currently have started at CERN. With the de- 
signed center-of-mass energy 14 TeV and luminosity L = 10 u cm~ 2 s~ 1 , the LHC gives 
access to high energy frontier at TeV scale and an opportunity to further improve the 
consistency test for the CKM matrix. The production rates for heavy quark flavours will 
be large at the LHC, and the bb production cross section will be of the order 0.5 mb, pro- 
viding as many as 0.5 x 10 12 bottom events per year (26[. In particular, the LHCb detector 
is designed to exploit large number of 6-hadrons produced at the LHC in order to make 
precise studies on CP asymmetries and on rare decays in 6-hadron systems. The other 
two experiments, ATLAS and CMS, are optimized for discovering new physics and will 
complete most of their B physics program within the first few years [26j, |27|. Obviously, 
the LHC has a great advantage over the current experiments on 6-hadrons[28(. 

In the present work, we have predicted possible large CP violating asymmetries in 
decay channel of B® — > K°p°(u) — > K°7t + tt~ via the p — u mixing. At the LHC, the 
6-hadrons are produced from the pp collisions. The possible asymmetry between the 
numbers of the 6-hadrons, Hb, and those of their antiparticles, Hb, has been studied in 



the Lund string fragmentation model and the intrinsic heavy quark model [29l. |30| . It has 
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been shown that this asymmetry can only reach values of a few percent. In our following 
discussions, we will ignore this small asymmetry and give the numbers of Hb&b pairs 
needed for observing the CP violating asymmetries we have predicted. These numbers 
depend on both the magnitudes of the CP violating asymmetries and the branching ratios 
of heavy hadron decays which are model dependent. For one (three) standard deviation 
signature, the number of H^Hf, pairs we need is 



where BR is the branching ratio for if& — > fp°. 

For central value of CP asymmetry in Eq. (!7H|) . we present the numbers of B S B S 
pairs for observing the large CP violating asymmetries at LHC For the channel 
B° s K°p°(u) ->■ K°n + n- , the numbers of B S B S pairs are 3.8 x 10 6 (3.4 x 10 7 ) for la 
(3a) si gna ture. At the LHC the average B S B S production is about 10% out of 10 12 bb 
events 29 ] . From Fig.l, one can see CP violating asymmetries vary sharply at small 



energy range, and reach peak value at a/s = 0.782 GeV. Hence, it is very possible to 
observe the large CP violating asymmetries in small energy range of p° ~ u resonance 
at the peak values of CP violating asymmetries from the LHC experiment. For the 
experiments, it is possible to reconstruction 7r + , 7r~ and K° mesons when the invariant 
masses of 7r + 7r~ pairs are in the vicinity of the u resonance. Therefore, it is very possible 
to observe the large CP violating asymmetries in B® — > K°p°(uj) — > fT°7r + 7r~ at the LHC. 



VIII. SUMMARY AND CONCLUSIONS 

In this paper, we have studied CP violation in B® — > K°p°(uj) — > K°7r + 7T~ . It has been 
found that, by including p — u mixing, the CP violating asymmetries can be large when 
the invariant masses of 7r + 7r~ pairs are in the vicinity of the u resonance. For the decay 
B® — > K°p°(u) — > K°7t + 7t~, the maximum CP violation can reach 46%. Furthermore, 
taking p — u mixing into account, we have calculated the branching ratio of the decay. 
We have also presented the numbers of B S B S pairs required for observing the predicted 
CP violation in experiments at the LHC. We have found the channel is the likely channel 
in which the large CP violating asymmetries may be observed at LHC. We expect that 
our predictions will provide a useful guidance for future investigations and experiments. 

In our calculations there are some uncertainties. We have worked in the QCD factor- 
ization which is expected to be a reliable approach in the heavy-quark limit. In the QCD 
factorization scheme, a s (rrib) and some l/m^ (annihilation) corrections are included. In 
this framework, there is cancellation of the scale and renormalization scheme dependence 
between the Wilson coefficients and the hadronic matrix elements. However, the QCD 
factorization scheme suffers from endpoint singularities which are not well controlled. The 
CP violating asymmetry depends on the unknown parameters which are associated with 
such endpoint singularities. The CKM matrix elements also lead to some uncertainties in 
the CP violating asymmetry. Uncertainties also come from the weak form factors associ- 
ated with the hadronic matrix elements. This lead to uncertain CP violating asymmetries 
in the QCD factorization scheme. This needs further detailed investigations. 
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